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Abstract 
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Manning-Rosen potential, the 'standard' Hulthen potential and the generalized Morse potential 
are briefly studied as special cases. It is observed that our analytical results are the same with the 
ones obtained before. 
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I. INTRODUCTION 



In this letter, we study the bound state solutions of a diatomic molecular potential having 
the form 

— fir —2fir 

V (r,M = -V° T ^ Fl ; + V 1Jr ^ V2 , (1) 

which has been firstly presented by Sun to fit experimental data of some diatomic molecular 
systems [1]. Analytical study of the above potential could be interesting since it involves 
several potential forms (for example, q = gives the Morse potential, q > corresponds to 
the 'generalized' Hulthen potential, etc.) meaning that we can simply extend the solutions 
to the ones of these special cases. 

The above potential is one of the central potentials which are a powerful ground for 
experimental and theoretical computations in different areas of physics such as in high energy 
physics where they were used to describe hadrons as bound states [2], in atomic physics where 
some important subjects such as binding energy and inclusive momentum distributions are 
studied by using of central potentials [3], in theoretical molecular dynamics model to study 
the intramolecular and intermolecular interactions and atomic pair correlation functions 
[4]. Moreover, the central potentials have been used in an important quantum mechanical 
problem which is also related with quantum information theory, the Fisher uncertainty 
relation and applied to the hydrogen atom and isotropic harmonic oscillator [5] and also 
for some theoretical calculations within the information theory to study some statistical 
quantities such as the Boltzmann-Shannon entropy [6]. The construction of an algorithm 
could be an interesting problem where the aim is to solve the radial Schrodinger equation 
(SE) for a given central potential V(r) numerically [7]. 

To our knowledge, the potential under consideration has been studied within the super- 
symmetric quantum mechanics [8] and in terms of Green's function [9] in non-relativistic 
domain. We search the bound state spectrum and the wave functions of the above potential 
by using an approximation instead of the centrifugal term in the same domain. We find 
an analytical expression for the energy spectrum and obtain the normalization constant by 
using some properties of the hypergeometric functions. Throughout this work, we restrict 
ourself to the cases where q > 1 and q = and give our numerical results for two diatomic 
molecules for different values of quantum number pair (n,£). 
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II. ENERGY SPECTRUM AND WAVE FUNCTIONS 



The radial Schrodinger equation is written [10] 



d 2 R{r) 
dr 2 



+ 



2m rrn . £(£+1) 

w [E ni - V(r)} - 



R{r) = , 



(2) 



where £ is the angular momentum quantum number, m is the particle mass, V(r) is the 
central potential and E n £ is the non-relativistic energy. Inserting Eq. (1) into Eq. (2) gives 



d 2 R(r) 
dr 2 



+ 



2mV l 



1 



+ 



2mV n 



1 2mE ne £(£+1) 



R(r) = . 



(3) 



k 2 (e^-g) 2 ' ti z e? r -q ' %- r 

where Vb, Vi, /? and g are real parameters defined by V\ = A)(e M — q), V = 2 VI, /3 = /i/r , 
where _D is the depth of the potential, r is the equilibrium of the molecule and q is the 
shape parameter. 

We use the following approximation [11] instead of the centrifugal term among the others 
[12-15] to obtain an analytical solution of Eq. (3) 
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(4) 



V (e? r - q) 2 ' 

Defining a new variable z = ge _/3r and taking a trial function as R(z) = z Al (l — z) A2 4>(z) 
and with the help of Eq. (4), Eq. (3) turns into 



z(l - z)*¥£ + [1 + 2 Al - (1 + 2M + 2A 2 )z] ^ 



dz 2 



where we set the parameters 



dz 



+ 



-2A X A 2 -A\ + 



2 2mVi 2mVo 



<P(z) = 0, 



2mE, 



n 



By using the abbreviations 



c = 1 + 2A 1 , 



b = A 1 + A 2 + x lA 2 + ^- 2 (v 1 + ^ l , 



a = Al + A 2 - x lA 2 + ^(v 1 + ^ , , 



(5) 

(6a) 
(6b) 

(7a) 
(7b) 

(7c) 
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Eq. (5) becomes an equation having the form of the hypergeometric-type equation [16] 



z{\ - z)<t>"{z) + [c - (a + b + l)z\4>'(z) - ab(f)(z) = , 



(8) 



whose solution is the hypergeometric functions 



<p(z) ~ 2 F 1 (a,b;c;z) . 



(9) 



When either a or b equals to a negative integer —n, the hypergeometric function <p{z) can 
be reduced to a finite solution. This gives us a polynomial of degree n in Eq. (9) and the 
following quantum condition 



Ai + A 2 - a \A\ 



2m 

W 1 



v 1 + ^ 

Q 



-n, 



(10) 



which gives the energy values of the two-term potential for any ^-values 



2*2 



E„e — 



(3 2 h 
2m 



TV 



+ (2n+l)(A' 2 + i)+i 



2n + 1 + 2A' 



(11) 



where 



A' 



1 £(£+1) 2mV 1 



(12) 



By using Eq. (10) we obtain the total wave functions 

R(z) = Nz Al (l - z) M 2 F 1 {-n, n + 2A 1 + 2A 2 ; 1 + 2A X ; z) . 



(13) 



where iV is the normalization constant and will be derived in Appendix A. 

We summarize our numerical results in Table 1 and 2 where the computations are made 
for two diatomic molecules, namely H 2 and LiH . The values of potential parameters we used 
for these molecules are as follows [13]: D = 4.744600 eV, r = 0.741600 A, m = 0.503910 
amu, /i = 1.440558 and E = U 2 /(mrl) = 1.508343932 x 10~ 2 eV for H 2 molecule and 
D = 2.515287 eV, r = 1.595600 A, m = 0.8801221 amu, /i = 1.7998368 and E = 
1.865528199 x 10~ 3 eV for LiH molecule [17]. It is seen that the energy values decrease while 
the values of the quantum numbers increase and the energy eigenvalues are also inversely 
proportional with the shape parameter for each of the molecules. 

Now we intend briefly to study some special cases whose energy eigenvalue equation 
obtained from Eq. (11) by suitable choices of the potential parameters. 



1. Manning-Rosen Potential 



The Manning-Rosen potential can be written as [18] 

Ah 2 1 ~'~ ^ 2 



V(r) 



2mb 2 e r / b - 1 



I a(a — l)h 



2mb 2 



[e'-l h - l) 2 ' 



(14) 



If we write our parameters as Vo = V\ — a{ % h ^ \P> = \ and q — 1 then we obtain the 
energy eigenvalues of the Manning-Rosen potential 

+ (2n + 1) (j + ^I + ^+l) + a (a-l) ) + ^ + 1) - A 



Enl ~ ~2~b 2 



n 



2n + 1 + 2y / i + £(£+l) + a(«-l) 



(15) 



which is the same result obtained in Ref. [19]. The normalization constant in Eq. (13) is 
obtained from Eq. (A8) as 

1/2 



N 



g(A?\ A$\k)g{A™, A%>, I) 2 F 1 {-2A£>, 1 + 2A\ L> + k + l;2 + 2A\ L > + k + Z; 1) 



(1) 



(16) 



where = ^-2mE ni b 2 /h 2 and = (1/2) (\ + yj 1 + U{t + 1) + 4ma(« - l)/h 2 ^ . 

We summarize our numerical results obtained from Eq. (15) in Table 3 where we set the 
parameters as A = 2b and a = 0.75 to compare the results with the ones given in Ref. [19]. 
Please note that the parameter D Q used in Ref. [19] is zero in the present work since our 
approximation used for the centrifugal term is different from the one used in Ref. [19] where 
energy eigenvalues are computed in atomic units. 



2. Standard Hulthen Potential 

Eq. (1) gives the standard Hulthen potential for V± — and q — 1 

V{r) = -V 



1 - e-f* ' 

and we obtain the energy eigenvalues from Eq. (11) 



p 2 h 

' 2m 



(n + £)(n + £ + 2) + l- 2 ^ 



2(n + l + £) 



(17) 



(18) 
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and the normalization constant of the corresponding wave functions from Eq. (A8) 



N = 



T(2 + 44 2) + k + i) 



1/2 



(19) 



_g(A u A { 2 2 \k)g(A 1 , A?\i)(l + 2Af + k + £)T(1 - 2A 1 + 4^ 2) )_ 

where A^p = 1 + £ . Choosing the parameters as (3 = ^ and V = a gives the following 
expression (m = H = 1) 



1 



nl 



2a 2 



{n + £)(n + £ + 2) + 1 -2aa 2 



(20) 



2(n+ 1 +£) 

which is the same result given in Ref. [20]. The standard Hulthen potential in Eq. (16) 
could gives the Coulomb potential for (3r <^ 1 

Ze 2 



V(r) 



(21) 



where we set V = Ze 2 /3 . We obtain the energy spectrum of the Coulomb potential from 
Eq. (18) (m = h = e = l) 

Z 2 

Ene -~2(n + l+£) 2 - (22) 
which is the same result obtained in Ref. [19]. The normalization constant of the corre- 
sponding wave functions is given with the help of Eq. (19) under the above assumptions. 

The numerical energy values of the Hulthen potential obtained from Eq. (18) are placed 
in second part of Table 3 for different quantum number pair (n, £) in atomic units. We 
choose the parameters as V = (5 = 5 as in Ref [21]. 

3. Generalized Morse Potential 



We obtain the generalized Morse potential for the limit q — > in Eq. (1) 

V(r) = V ie - 2(3r - V Q e-P r , (23) 



which gives the following equation for s-waves 
dr 2 



d*_ _ 2mV 1 o r 2mV 2/3r 
lr ' Ir + 



2mE n A 

— ) R ( r ) = > 



ft / 



(24) 



Defining a new variable z = e /3r and taking the wave function of the form R(z) 
e~ Blz / 2 z B2 / 2 (f)(z), we obtain 



d 2 (f){z) . . d(j)(z) 

z—^ 1 + (1 + B 2 - B lZ ) + 



dz 2 



dz 



B L B 2 B l 2mV 
~2 2 + (3^h 2 



<l>(z) = 0, 



(25) 
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where B\ = and B\ = Using a new variable y = B 1 z gives 



d My) , „ x d6(z) 



dy 2 dz 
which is the Laguerre differential equation 



B 2 1 2mV 



2 2 B^h 2 



</%) = 0, (26) 



xy" + (a + l- y)y' + ny = . (27) 

where the factor n should be zero or a positive integer to get a polynomial solution [22]. So, 
the solution of Eq. (25) are given in terms of the Laguerre polynomials as 

Hv) ~ Ll{y) , (28) 

where a = B 2 and n = — ^ — | + • We get the total eigenf unctions of the Morse 

potential 

R(z) = Ne- B ' z/2 z B2/2 Ll(B lZ ) . (29) 

and the energy eigenvalues 

£ --^r +i -^| • (3o) 

We present the numerical energy values of the Morse potential obtained from the above 
equation in Table 3. We give the results for H 2 molecule (in eV) by taking the same param- 
eter values to obtain the results given in Table 1 and by setting the potential parameters as 
V = 2D , V, = D . 

Using the following representation of the Laguerre polynomials [22] 

m*) = ±(-v( n+v )i. (3D 



k=o \n — k 

the normalization condition is written as 



\N\ 2 g(n)g{m) [ e BlZ z B * +2k dz = 1 , 
Jo 



(32) 



where 



k=0 \ It- — rv / k=0 



Changing the variable t = B\z in Eq. (32) gives 

\N\ 2 g(n)g(m)B^ 1+B2+2k) f e'H^dt = 1 , (34) 

Jo 

which includes the incomplete Gamma function defined as [22] 

r i 

7(0,0;)=/ t a - 1 e- t dt = -x a e- x 1 F 1 (l;l + a;x), (35) 
Jo « 

Finally the normalization constant is obtained as 



N 



g(n)g(m) 1 F 1 (l;l + n;l) 



, tt = l + B 2 + 2k. (36) 



III. CONCLUSIONS 



We have studied the approximate bound state solutions of the radial SE equation for 
a two-term potential. We have obtained the energy eigenvalues and the corresponding 
normalized wave functions approximately in terms of the hypergeometric functions. We 
have presented our numerical results of the energy eigenvalues of two diatomic molecules in 
Tables 1 and 2. We have also studied the analytical bound state solutions of the Manning- 
Rosen potential, the 'standard' Hulthen potential and the generalized Morse potential as 
special cases. We have observed that our all analytical results are the same with the ones 
obtained in the literature. We have also summarized some numerical results of the energy 
eigenvalues of the above three potentials in Table 3 and observed that our results are good 
agreement with the ones obtained before. 
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Appendix A: Normalization Constant 

The wave functions in Eq. (13) is 

R{z) = Nz Al (l - z) M 2 F 1 (-n, n + 2A l + 2A 2 ; l + 2A 1 ;z), (Al) 
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which is written in terms of the new variable z = g£ (0 < £ < 1) 

R(qO = Nq A ^ Al (l - qi) M tF^-n, n + 2A 1 + 2A 2 ; 1 + 2A ±] q£) , 
The normalization condition f 1 \R(q£)\ 2 d£ — 1 gives 

|Af q 1+2Al f e Al (l - q0 2M hF^-n, n + 2A 1 + 2A 2 ; 1 + 2A X - g£)f dx = 1 
Jo 

Using the representation of the hypergeometric functions [22] 

tp < u \ (- n )k{b)k k 

Eq. (A3) becomes 

|Af q 1+2M g(A u A 2 , k)g{A u A 2 , 1) f ^ +k+l (l - qtf^dt = 1 , 

Jo 

where (_ n ) fc = (_i)*( n - + l) fc = (-l) fc and 

oM , m f- (-n) fc (n + 2A 1 + 2A 2 ) fc fc 

and g(A h A 2 , 1) = g(A u A 2 , k ->■ Z) . 

Using the following identity for the hypergeometric functions [22] 

2*1 *) = r(^)r(y - jf " ^"'(i " > 

we obtain the normalization constant from Eq. (A5) 



N 



1 



g i+2A^ (Ai; A ^ k )g(A u A 2 , 1) 2 Fi(-2A 2 , 1 + 2Ai + fc + Z; 2 + 2A 2 + k + Z; g). 
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TABLE I: Energy eigenvalues of the J 


H 2 molecule in eV (E n g < 0). 




n t 


q= 1.25 


q = 1.50 


q = 1.75 





3.8099200 


2.3409000 


1.5021100 


1 


2.6465800 


1.4911800 


0.8623870 


1 


2.6295900 


1.4804000 


0.8551840 


2 


1.7526700 


0.8721010 


0.4265830 


1 


1.7394900 


0.8642480 


0.4217640 


2 


1.7133000 


0.8486670 


0.4122170 


3 


1.0823900 


0.4428370 


0.1570130 


1 


1.0724800 


0.4374880 


0.1542230 


2 


1.0528200 


0.4269000 


0.1487240 


3 


1.0237200 


0.4112900 


0.1406770 


4 


0.5991580 


0.1711030 


0.0242272 


1 


0.5920830 


0.1679170 


0.0231832 


2 


0.5780760 


0.1616440 


0.0211659 


3 


0.5574210 


0.1524740 


0.0183168 


4 


0.5305400 


0.1406920 


0.0148458 


5 


0.2734990 


0.0311495 


0.0049878 


1 


0.2689040 


0.0298494 


0.0054591 


2 


0.2598410 


0.0273349 


0.0064646 


3 


0.2465630 


0.0237765 


0.0081293 


4 


0.2294460 


0.0194276 


0.0106393 


5 


0.2089850 


0.0146225 


0.0142399 
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TABLE II: Energy eigenvalues 


of the LiH molecule in eV (E ni < 0). 




n t 


q= 1.25 


q = 1.50 


q = 1.75 





3.5677900 


2.3156300 


1.5836900 


1 


3.0290300 


1.9055500 


1.2595300 


1 


3.0250700 


1.9029400 


1.2577000 


2 


2.5506900 


1.5476900 


0.9820380 


1 


2.5471400 


1.5453900 


0.9804600 


2 


2.5400700 


1.5408100 


0.9773080 


3 


2.1274300 


1.2372900 


0.7468470 


1 


2.1242700 


1.2352900 


0.7455040 


2 


2.1179600 


1.2312900 


0.7428200 


3 


2.1085100 


1.2253000 


0.7388050 


4 


1.7544800 


0.9701380 


0.5500770 


1 


1.7516800 


0.9684030 


0.5489500 


2 


1.7460700 


0.9649380 


0.5466990 


3 


1.7376900 


0.9597550 


0.5433330 


4 


1.7265600 


0.9528690 


0.5388640 


5 


1.4275700 


0.7424420 


0.3882750 


1 


1.4251000 


0.7409570 


0.3873480 


2 


1.4201500 


0.7379920 


0.3854990 


3 


1.4127600 


0.7335570 


0.3827350 


4 


1.4029300 


0.7276680 


0.3790660 


5 


1.3907100 


0.7203470 


0.3745090 
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TABLE III: Energy 


eigenvalues obtained from Eq. 


(15), Eq. (18) and Eq. 


(30). 


Manning-Rosen potential 


n 


I 


1/6 




Present Work 


Ref. [19] 


2 


1 


0.025 




-0.1205793 


-0.1205279 






0.050 




-0.1084228 


-0.1082170 






0.075 




-0.0969120 


-0.0964490 






0.100 




-0.0860470 


-0.0852240 


Hulthen potential 


n 


I 


5 




Present Work 


Ref. [21] 





1 


0.025 




-0.1128130 


-0.1127600 






0.050 




-0.1012500 


-0.1010420 






0.075 




-0.0903120 


-0.0898450 






0.100 




-0.0800000 


-0.0791700 






0.150 




-0.0612500 


-0.0594950 


Morse potential 


n 








Present Work 


Ref. [17] 











-4.476013 


-4.476013 


1 








-3.962315 


-3.962315 


2 








-3.479919 


-3.479918 


3 








-3.028824 


-3.028823 


4 








-2.609030 


-2.609029 


5 








-2.220537 


-2.220536 
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